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Abstract. A consistent analytical approach for calculation of the 
quasiclassical radial dipole matrix elements in the momentum and 
coordinate representations is presented. Very simple but relatively 
precise expressions for the matrix elements are derived in both rep- 
resentations. All analytical expressions contain only one special 
function - the Anger function and its derivative. They generalise 
and increase the accuracy of some known quasiclassical expressions. 
The small difference between the two forms of the expressions for 
the dipole matrix elements indicates to the applicability of the sim- 
ple expressions given by the consistent quasiclassical approach even 
for low atomic states. 
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1. Introduction 

Calculation of transition probabilities, oscillator strengths and dipole matrix 
elements for the atomic transitions is stimulated of investigations in spec- 
troscopy, plasma physics, chaotic dynamics of the nonlinear systems and 
other fundamental and applied fields. The accurate calculations of these 
atomic characteristics are very time consuming and in some cases involve 
certain difficulties, especially for transitions between high states. Therefore, 
recently new attention has been paid to quasiclassical calculations of the 
dipole matrix elements and oscillator strengths for atomic nl — > n'l' transi- 
tions (Heim et al 1989, Delone et al 1989, 1994, Kaulakys 1991, D'yachkov 
and Pankratov 1991, 1994, Pankratov and Meyer-ter Vehn 1992 and Nana 
et al 1995). In the review paper by Delone et al (1994) the set of analyti- 
cal formulae for the dipole matrix elements between quasiclassical states is 
presented. However, this review paper does not reflect the results of the 
last five-year period by Heim et al (1989), Kaulakys (1991), D'yachkov and 
Pankratov (1991, 1994) Pankratov and Meyer-ter Vehn (1992) and others 
and contains some inconsistencies. As a result the main formulae of the 
paper contain terms with erroneous sign and do not agree with the exact 
results, e.g. when s = v' — v is an integer number and when s — ■> with v 
and v' being the effective principal quantum number of the initial and final 
state, respectively. 

Here we present a consistent analytical approach for the quasiclassical 
dipole coupling of the electron with the electromagnetic field, taking into 
account peculiarity of the radial quasiclassical matrix elements: radial wave 
functions of the initial and final states for transition with Al ^ correspond 
to the different effective Hamiltonians and, therefore, we can not use the 
usual correspondence between the quasiclassical matrix elements and the 
Fourier components of the appropriate classical variable (see e.g. Landau 
and Lifshitz 1985). So we should start from the definition of the quantum 
matrix elements using the quasiclassical radial wave functions. 

The direct way of coupling a radiation field to the electron Hamiltonian is 
through the A • p interaction (see, e.g. Bethe and Salpeter 1957 and Landau 
and Lifshitz 1985) where A is the vector potential of the electromagnetic 
field and p = —ihV is the momentum operator. When the radiation wave- 
length is long compared with atomic dimension, as is the case for optical or 
microwave transitions, the variation of the vector potential within the atom 
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can be neglected. In this case the electric field is also uniform over the extent 
of the atom. So that the electric dipole moment is the only atomic multipole 
coupled to the field in the multipolar approximation to the interaction oper- 
ator and the interaction of the electron with the field may also be expressed 
through the E • r term, where E is the electric field strength and r is the 
electron's coordinate. As a matter of fact, the two forms of the interaction 
Hamiltonian in the dipole approximation are equivalent due to the gauge 
invariance of the field. 

Usually one calculates the radial dipole matrix elements in the coordinate 
form (see Heim et al 1989, Delone et al 1989, 1994, Kaulakys 1991, D'yachkov 
and Pankratov 1991, 1994, Pankratov and Meyer-ter Vehn 1992, Nana et al 
1995 and references therein). For the precise wave functions of the | a) and 
| b) states the relation between the matrix elements of p and r operators 

Pa,b = -mwj/^ (1) 

holds. Here m and Ub A = (E b — E a )/fi are the electron mass and the angular 
frequency of the corresponding transition, respectively. However, when one 
uses approximate wave functions in the evaluation of the matrix elements, 
the length and velocity forms in general yield different results and the re- 
lationship (1) holds only roughly or even fails (see, e.g. Crossley 1969). If 
asymptotically correct wave functions are used, the r-form of the dipole ma- 
trix elements is preferable as it stresses the contributions to the integral from 
the large r region. On the other hand, for the wave functions correct at small 
and medium r the velocity form should be used as it puts more weight on 
the integral at small and medium r. The quasiclassical wave functions in the 
classically allowed region of motion are not correct asymptotically as well as 
for very small r but are relatively accurate for the medium r between the two 
turning points of the classical orbit. Therefore, it is likely that the velocity 
form of the quasiclassical matrix elements is as much (or, maybe, even more) 
accurate as the length form. 

2. Radial dipole matrix elements for any spherical potential 

Calculation of the angular part of the matrix elements is a simple problem 
and, therefore, we restrict ourselves to the radial part. The radial dipole 
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matrix element in the momentum representation is given by Q 
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were P n i{r) = rR n i(r) is the solution of the radial Schrodinger equation, 
uj = E n iy — E n i is the transition frequency and Z max = max(/,/ ± 1). The 
quasiclassical radial wave function P nl can be expressed as 
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in the classically allowed region of motion and some exponentially decreasing 
function outside this region. Here T is the period of classical rotation, the 
radial velocity of the electron v r is given by 



v r (r) 

and the phase $ n / is defined as 
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with U (r) being the effective potential which defines the motion of the Ry- 
dberg electron and roots of the radial velocity v r (r), r\ and r2, are the two 
classical turning points given by v r {r\) = v (r 2 ) = 0. 

Further calculations of the matrix elements (2) are similar to those by 
Heim et al (1989) and Kaulakys (1991). Substituting equation (3) into equa- 
tion (2) and neglecting the rapidly oscillating sine and cosine of the sum of 
the phases $ n z and <& n 'i' one finally finds that 
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Further in the paper we will use the atomic units, h — m — e = I. 
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where ip is the polar angle while x and y are Cartesian coordinates of the 
electron, points denote the derivatives with respect to time and T c is some 
mean period to be defined later. In the derivation of equation (6) we have 
used the fact that r" 1 /^ = r~\l + \) = r~ l {l + V + l)/2 = r _1 (rV) = rtp 
and that, according to equations (4) and (5), linear expansion of the difference 
of the phases in powers of u and Al = I' — I is 

A$ = $ n /,,-$ n ,~wf-AZp + .... (7) 

Integration of equation (6) by parts yields to the r-form of the radial 
dipole matrix element 



T c /2 

Rtl ±l = {nl | r | n'l ± 1) = ^= / r(t) cos A*(t)]dt 



o 

T c /2 



2 f 

= ^Jff} J [x(t)cosut±y(t)sinut]dt. (8) 

Here we have used the quasiclassical quantisation conditions (see Kaulakys 
1991 for details) and facts, that 

A$(0) ~ 0, A$(T c /2) = An r 7r, uoT c /2 = Ann, 

r (0) = x(0) = n, y(0)=y(T c /2) = 0, <p(0) = 0, V (T c /2) = tt (9) 

where n r = n — I — 1 is the radial quantum number and An and An r are 
integers. Note, that equation (8) may also be derived from the definition of 
the radial dipole matrix element in the coordinate representation and using 
expression (3) for the quasiclassical radial wave function (see Kaulakys 1991). 

It follows from equations (6)-(8) that only for Al = the quasiclassical 
matrix elements of the module of the radius vector r coincide with the Fourier 
components of the classical variable r(t) and may be expressed through the 
Fourier components of the classical variable f(t), i.e. 

T c /2 

2 f If 

(nl\r\ n'l) = -j== J r(t) cos utdt = -j== j r(t)e' tujt dt 
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For the dipole transitions with Al = ±1 radial wave functions P n \ and 
P n >i> of the initial and final states are solutions of the radial Schrodinger 
equation with different effective potentials V e //(r) —U(r) + (l+ |) 2 /2r 2 and 
Vejf(r) = U(r) + (V + |) 2 /2r 2 , respectively. This results to the additional 
term in the phase difference (7) and second term in equations (6) and (8). 

Note, that the radial dipole matrix elements D^ l±1 and R^ l±1 according 
to equations (6) and (8) may also be expressed as 



It is of interest to note the connection between the expressions for the 
dipole matrix elements (6) - (10) and the energy change of the classical 
atom in a monochromatic field (see Gontis and Kaulakys 1987 and Kaulakys 
1991). The mapping equations of motion for the classical hydrogen atom 
in an oscillating electric field derived on the bases of the velocity form of 
interaction are suitable for investigation of transition to chaotic behaviour 
and ionisation of Rydberg atoms even in the low frequency field when the 
strength of the external field is comparable with the Coulomb field of the 
atom (Kaulakys and Vilutis 1995). 

Until now we did not use the explicit form of the potential U(r). There- 
fore, equations (1)— (12) are valid for any spherical potential. Further we will 
apply this theory for the non-hydrogenic atoms. 

3. The non-hydrogenic atom 

The potential U (r) which defines the motion of the Rydberg electron of the 
non-hydrogenic atom or ion may be expressed as a sum of two terms - 
the Coulomb potential —Z/r of the ion core with charge Z = Z ion + 1 and 
the perturbation potential AU(r) due to the deviation from the Coulomb 
approximation, i.e. 




(11) 



and 




(12) 



U(r) = -Z/r + AU(r). (13) 

For the non-hydrogenic atom the perturbation potential AU (r) is short-range 
and results to the non-Coulomb scattering phase Si related with the quantum 
defect pLi by the relationship Si = 7r/z;.The energy E nl of the | n/)-state is 
related to the effective principal quantum number v and the quantum defect 
Hi in the usual way 

E nl = -Z 2 /2v 2 , v = n-m. (14) 

Significant contribution to the integrals (6) and (8) arise from the regions 
with relatively large r where the potential U (r) is well represented by the 
Coulomb potential V (r) = —Z/r. The additional potential AU(r) results 
mainly to the non-Coulomb scattering phaseshift Si (see Kaulakys 1991). 
Thus, the phase (5) and the phase difference (7) in the region of the main 
contribution to the dipole matrix element may be represented as 

$ni(r) = f v c r .{r)dr + S t - ]ir (15) 

Jr\ 4 

and 

A$ ~ out - Alip ~ A + uot c - Alip c (16) 

where v c r is the radial velocity for the Coulomb potential, A = Sy —Si and it is 
convenient to introduce the parametric equations of motion for the Coulomb 
potential 

r = (u 2 jz) (1 - ecosO , t c = (yljZ 2 ) (f - esinf) 

x=(ul/z){co^-e), y=(^)(l-e 2 ) 1/2 sin^ (17) 

cos £ — e \ 
1 -ecosU ' 
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Here e denotes the eccentricity of the classical orbit with the period T c = 
27rz/^/Z 2 and the turning points r\ 2 = (u 2 /Z) (1 q= e) . The centered effective 
principal quantum number v c is defined from the requirement that the phase 
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difference (16) at the turning point r2 has to be in consistence with the 
quasiclassical quantisation conditions, i.e. 



A$(r 2 c )=A + -a/T c 



Aln = An r 7i, 



(18) 



which, together with the relationship An r 7r = Avtx + A — Aln, results to the 
expression (see also D'yachkov and Pankratov 1991, 1994) 



Z 2 Av 2 (vv'\ 



v + v' 



(19) 



More precisely, the non-Coulomb phase shifts Si and Sy and, consequently 
the phase difference A = 8y — Si are some functions of the r and time t (see 
Kaulakys 1991). However, for the non-hydrogenic atoms the phase difference 
A (r) increases in the region r ~ r\ ~ r\ from A = to A = Sy — Si — const 
in the very narrow interval of the coordinate r, while the main contributions 
to the radial integrals (6) and (8) occur at large distances, r ~ v 2 . Thus, in 
equations (6) and (8) the phase out may be replaced by the expression A+ut c , 
while ip ~ ip c . Keeping this in mind and substituting equations (16)-(19) into 
equations (6) and (8) we have 
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sin 7rs, s = Av = v' — v. 



Here J_ s (z) and J'_ s (z) are the Anger function defined as 



1 r w 

— / cos (s^ + z sin ^) 
% Jo 



(20) 

(21) 
(22) 
(23) 

(24) 



and its derivative with respect to the argument z, respectively. Note to the 
properties of the Anger function: J_ s (z) = J s (—z) and J_ s (z) = —J s (—z) 
which result to the symmetry of the matrix elements (21) and (23) 
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D ± (e,-s) = D T (e,s). (23') 

To the best of our knowledge equations (20) and (21) are derived for 
the first time while equations (22) and (23) within the accuracy of fac- 

tor [y c j V vv' j coincide with the equation (16) in the paper by Kaulakys 
(1991) and are close to the corresponding expressions given by D'yachkov 
and Pankratov (1991) and by Pankratov and Meyer-ter-Vehn (1992) (see also 
Nana et al 1995). Note, that for the first time the dipole matrix elements in 
the coordinate representation have been expressed through the Anger func- 
tion and its derivative (however, with some erroneous signs) by Davydkin 
and Zon (1981). 

In the derivation of equations (22) and (23) one integrates equation (8) 
by parts using the approximate expression ut ~ A + ut c . This yields to 
equation (6) with the additional term 

RT 1 - Dt! ±l = ^L= sin A = (-l) An ^-tf^ (25) 

nl nl ujVTV V ; Z (z/z/') ITS 

Thus, the difference between the quasiclassical dipole matrix elements in the 
r-and ]?-forms, the second term in equation (23), results from the replacement 
of the electron's motion in the effective potential U (r) by the motion in the 
Coulomb potential with the additional phase 8\. Therefore, the quasiclassical 
radial dipole matrix element in the coordinate representation (22)-(23) may 
contain some additional inaccuracy. On the other hand, the difference (25) 
between two forms of the dipole matrix elements may be as a criterion of the 
exactness of the quasiclassical approximation. As a rule, the additional term 
(25) is small because of the small factor (1 — e) for states with low / and 
of the small factor sin7rs for states with larger / but small quantum defects 
and, consequently, with s = Au ~ An close to the integer. In expansion of 
the dipole matrix elements in terms of a = y/l — e 2 = (/ + /' + 1) /2u c (see 
Kaulakys 1991 for analogy) 



p s 



J'_ s (s) ± a [3^ s (s) - — — J + — (^r s (s) + 



2 sin 7ts 



TVS 



(26) 



a 2 sin tt.s 



D * = D ' + y— < 27) 
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this additional term makes up less than half of the third, proportional to 
the ct 2 , expansion term. This indicates to the relatively high (up to some 
percents) accuracy of the very simple quasiclassical approximation (20)-(23) 
for the dipole matrix elements even for the low atomic states. The extensive 
analysis of such approach and comparison with the numerical Hartree-Fock 
calculations will be presented elsewhere. Here we will present only the lim- 
iting forms of the dipole matrix elements. 

4. Special cases of the parameters 

Using the expansions of the functions J_ s (es) and J'_ s (es) in powers of s 
(Kaulakys 1991) 

J s (es) ~ 1 - Qtt 2 + e + ^e 2 ) s 2 

J ' s (es) ~ - (l + ie) s, s « 1 (28) 
we have from equations (20)-(23) 

D ni i±i = (-i) An+1 y i 1 + r) ^ 

Rtr = {-l) An+1 ~e s«l. (30) 
For hydrogenic atom with n' — n equation (30) results to the exact expression 

RT 1 = -~e. (31) 

Substitution of the asymptotic, s >> 1, forms of functions J_ s (s) and J'_ s (s) 
(Kaulakys 1991) 

2a / 1 \ 

J ~ s (s) = Vs^~ 3 c ° s V s ~ 6 V a ^ °' 447 
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into equations (26) and (27) yields 
26 / 1 
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(33) 



From equation (33) the Bethe rule can be seen: principal and orbital quantum 
numbers change prevailing in the same direction but only when a and / are 
not small and cos (ns + |7r) and cos (its — |7r) are of the same sign, e.g. 
when s is close to the integer. In the later case (s = An) we have the 
improved result of Goreslavsky et al (1982) for the removed states 
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On the other hand, for large s = An the Anger function and its derivative 
may be expressed through the Airy function and its derivative or through 
the McDonald functions. As a result we have from equations (20)-(23) 
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and K v (z) are the McDonald functions. 

For 1 — e 2 = a 2 1 it yields from equations (35) and (36) 
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In the limit ^a 3 An 1 equation (37) results to the expression (34). 

The dipole matrix elements for transitions between states with the large, 
/ ~ n, orbital quantum numbers, as follows from equations (20)-(23) or (35) 
and (36) when e — > 0, are exponentially small. Moreover, the Bethe rule 
in this case is enhanced: the transitions with the change of principal and 
orbital quantum numbers in the opposite directions are strongly suppressed 
in comparison with transitions, when n and / change in the same direction. 
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Thus, the very simple expressions (20)-(23) cover all known quasiclassi- 
cal non-relativistic results for the dipole matrix elements. They generalize 
and increase the accuracy of some earlier derived expressions. Extension of 
the present approach to the continuum states is rather straightforward (see 
Kaulakys 1991 for analysis in the r-representation). 

5. Conclusions 

Consistent analytical approach for calculation of the quasiclassical radial 
dipole matrix elements in the momentum and coordinate representations is 
presented and very simple but relatively precise expressions for the matrix 
elements are derived in both representations. All analytical expressions for 
the quasiclassical radial matrix elements in both representations contain only 
one special function - the Anger function and its derivative. They general- 
ize and increase the accuracy of some known quasiclassical expressions. The 
small difference between the two forms of the expressions for the dipole ma- 
trix elements indicates to the applicability of the simple expressions given by 
the consistent quasiclassical approach even for low atomic states. 

It is important to note that the dipole matrix elements as the analytical 
functions (even for the hydrogenic atom) are expressed through the Anger but 
not through the Bessel functions. It is another thing that the Anger functions 
Jj, (z) of the integer order v — m coincide with the Bessel functions J m (z), i.e. 
J m (z) = (— l) m J_ m (z) = J m (z). Expression of the dipole matrix elements 
through the Bessel function J s (es) or through the Anger function of the 
positive order and positive argument J s (es) (see, e.g. Delone et al 1994) 
results to the erroneous limit when s — > and to another inaccuracies. 
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